RD-A141  522  ITERATIVE  COUPLING  OF  PARTITIONED  CLUSTERS  IN  THE  1/ 

SCATTERED-HAVE  HOLECULA  .  OJ)  MASSACHUSETTS  INST  OF  TECH 
CAMBRIDGE  CENTER  FOR  MATERIALS  SCI.  .  F  A  LEON  ET  RL. 
UNCLASSIFIED  15  MAV  84  TR-?  N00014-81-K-0499  F/G  7/4  NL 


OIIC  FILE  COPY 


AD- A 141  522 


a ArP  READ  INSTRUCTIONS 

rMVjC _ BEFORE  COMPLETING  FORM 

i2,  GOVT  ACCESSION  NO,  3-  REORI En  T's  c  at  Al_o"c  NUMBER 


14.  T|TLE  (and  Subtttlt) 


Iterative  Coupling  Of  Partitioned  Clusters 
In  The  Scattered-Wave  Molecular-Orbital  Method 


5.  Type  OF  REPORT  A  PERtOO  CO VEREO 


Interim 


S.  PERFORMING  ORG.  REPORT  NUMBER 


I?,  author^; 


F.  A.  Leon  and  K.  H.  Johnson 


S.  CONTRACT  OR  grant  NUMBER/*) 


N00014-/81-K-0499 


I*.  PERFORMING  ORGANIZATION  NAME  ANO  AOORESS 


Center  for  Materials  Science  and  Engineering, 
M.I.T.,  Cambridge,  Massachusetts  02139 


II.  CONTROLLING  OFFICE  NAME  ANO  AOORESS 


Office  of  Naval  Research 
Department  of  the  Navy 
Arlington,  Virginia  2??i 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  ft  WORK  UNIT  NUMBERS 

Task  No.  Nr  056-757 


12.  REPORT  DATE 

May  15,  1984 

13.  NUMBER  OF  PAGES 

18 


MONITORING  AGENCY  NAME  ft  ADDRESS/!/  dittarant  tram  Cantrelllnf  Ottlea)  )  IS.  SECURITY  CLASS,  (at  Ui la  raport) 


l«.  DISTRIBUTION  STATEMENT  (at  Oita  Raparl) 

Approved  for  publfc  release; 
distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  (at  Ilia  abstract  antarad  In  alack  20.  II  dlllarant  i 


ft*  SKCLASSI  FI  CATION/ DOWNGRADING 
SCHEDULE 


I  Rapart) 


DTIC 

[EEECTE] 

MAY  2  3  1984 


I  I  ft.  supplementary  notes 


>t.  ICIY  monO%  f'Cam/iwt  on  rnnfi<  tide  it  »ncenf  mtd  l4m*Htr  bf  block  wterj 

iterative  coupling 
/  partitioned  clusters 

(  scattered-wave  molecular-orbital 

20.  \ABSTRACT  (Canttnaa  an  raaaraa  MW  H  nagatOmr  and  tdantttr  fty  fttoaft  mm*m) 

^)The  Iterative  coupling  technique  Is  described  and  applied  to  the  parti¬ 
tioned  self-consistent-field  scattered-wave  method.  This  technique 
allows  the  electronic  coupling  of  a  central  cluster  of  atoms  to  Its 
environment  to  be  calculated  iteratively,  permitting  scattered-wave  molecular- 
orbital  calcluatlons  to  be  performed  on  larger  clusters  than  was  previously 
possible.  A  test  calculation  on  methanol  Is  reported, 


DO,j2Tts  1473  EDITION  OF  I  NOV  «S  IS  OBSOLETE  UnClaSSlf1@d 

s/w  0102-LP-oi  j-6401  unci  ass  iTi  ea 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  fWiwi  Oata  Mnimad) 


OFFICE  OF  NAVAL  RESEARCH 


Contract  N00014-81-K-0499 
Task  No.  Nr  056-757 


TECHNICAL  REPORT  NO.  7 


ITERATIVE  COUPLING  OF  PARTITIONED  CLUSTERS 
IN  THE  SCATTERED-WAVE  MOLECULAR-ORBITAL  METHOD 


by 

F.  A.  Leon  and  K.  H.  Johnson 


Center  for  Materials  Science  and  Engineering 
Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts  02139 


May  15,  1984 


Reproduction  In  whole  or  in  part  is  permitted  for 
any  purpose  of  the  United  States  Government 


"v 

■V 

>> 

m 

s 

,V  * 

■Si 

fcts 


Approved  for  Public  Release:  Distribution  Unlimited 


84  0  5  2l  22  9 


V  v 


^ V"  W ;% 


ITERATIVE  COUPLING  OF  PARTITIONED  CLUSTERS 
IN  THE  SCATTERED-WAVE  MOLECULAR-ORBITAL  METHOD* 


F.  A.  Leon  and  K.  H„  Johnson 
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Cambridge,  Massachusetts  02139 


Abstract 


The  iterative  coupling  technique  is  described  and  applied  to 
the  partitioned  self-consistent- field  Xa  scattered-wave  method. 

This  technique  allows  the  electronic  coupling  of  a  central  cluster 
of  atoms  to  its  environment  to  be  calculated  iteratively,  permitting 
scattered-wave  molecular-orbital  calculations  to  be  performed  on 
larger  clusters  than  was  previously  possible.  A  test  calculation 
on  methanol  is  reported. 
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I.  Introduction 

The  multiple-scattering  or  scattered-wave  method  [1],  in  conjunction 
with  a  density- functional  approximation  to  exchange-correlation  such  as 
the  Xa  approximation,  is  a  widely  used  tool  for  molecular-orbital  cal¬ 
culations.  This  is  largely  due  to  its  computational  efficiency  and 
facility  in  dealing  with  large  polyatomic  molecules,  especially  those 
containing  heavy  atoms  (e.g.,  transition-metal  complexes).  Nevertheless, 
there  are  a  wide  variety  of  interesting  electronic  structure  problems 
that  involve  a  number  of  atoms  greater  than  that  which  is  practical  or 
possible  with  the  standard  scattered-wave  method. 

Often,  the  constituents  of  these  large  clusters  or  molecules  can  be 
divided  into  two  groups  of  atoms  (i.e.,  two  subclusters) :  one  having 
the  most  scientific  interest  (the  "principal"  subcluster),  and  the  other 
(the  "environmental"  subcluster)  having  less  interest  except  in  modifying 
the  properties  of  the  first  group.  For  example,  a  biological  molecule 
may  have  an  active  enzymatic  site  consisting  of  a  transition  metal  and 
Its  nearest  neighbors.  The  enzymatic  properties  of  this  site  might  be 
modified  by  other  subgroups  which  are  physically  more  distant.  In  this 
case  It  would  be  desirable  to  be  able  to  approximate  the  electronic 
coupling  between  the  more  distance  subgroups  and  the  complex  consisting 
of  the  transition  metal  and  its  first  nearest  neighbors.  For  a  large 
molecule,  such  an  approximation  would  allow  calculations  to  be  performed 
which  might  otherwise  be  Impossible.  Another  example  would  be  the  problem 
of  calculating  the  electronic  structure  of  a  cluster  that  represents  a  part 
of  a  solid  or  interface,  In  which  the  effects  of  "embedding"  the  cluster 
are  to  be  Included. 


In  this  paper,  a  method  for  performing  this  type  of  approximate 
calculation  is  described,  and  results  of  numerical  calculations  are 
reported.  The  basis  for  the  calculations  is  the  partitioned  form  of 
the  scattered-wave  method,  developed  by  Kjellander  [2-4],  in  which  the 
atoms  of  the  molecule  or  cluster  are  distributed  among  spherical  sub¬ 
clusters.  The  full  secular  matrix  of  the  partitioned  problem  is  then 
contracted,  leading  to  a  much  simpler  secular  problem.  The  molecular 
orbitals  of  the  approximately  coupled  scheme  are  found  by  an  iterative 
scheme  which  rapidly  converges  to  the  "exact"  molecular -orbital  energies 
and  wavefunctions;  this  result  is  demonstrated  by  a  test  calculation 
using  the  methanol  molecule. 

II.  Method 

In  the  partitioned  form  of  the  scattered-wave  method,  the  atoms  of 
the  molecule  or  cluster  are  divided  into  subclusters,  and  each  subcluster 
is  surrounded  by  a  partitioning  sphere.  While  in  the  standard  scattered- 
wave  method  the  interatomic  charge  and  potential  are  averaged  over  the 
entire  Interatomic  region,  in  the  partitioned  case  this  averaging  is 
performed  separately  for  each  subcluster.  Thus,  in  the  partitioned 
method,  the  constant  interatomc  potential  for  each  subcluster  is 
different,  leading  to  a  more  realistic  model  potential. 

Consider  a  partitioned  scattered-wave  calculation  like  the  one 
discussed  In  the  Introduction,  in  which  the  cluster  can  be  divided 
into  a  principal  and  an  environmental  subcluster.  The  secular 
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equation  of  this  problem  could  be  written  in  the  following  schematic 
form: 


In  this  equation,  the  matrix  T  represents  the  scattering  matrix  for  the 

a 

principal  subcluster  (including  scattering  from  the  sphere  surrounding  it), 
and  Tb”^  represents  the  scattering  matrix  for  the  environmental  subcluster; 
Gab  and  Gbfl  are  matrices  which  represent  the  Green's  function  that  describes 
the  propagation  of  the  electrons  between  the  two  subclusters;  and  Ag  and 
are  column  vectors  containing  the  atomic  and  partitioning  sphere  coefficients 
of  the  wavefunction.  (The  detailed  definitions  of  these  matrices  and  column 
vectors  will  be  given  in  a  future  paper,  or  can  be  found  by  studying 
Refs.  [11  and  [21.)  Note  that  (1)  Is  a  linear  system  for  the  coefficients 
Aa  and  Afa;  solutions  are  found  by  varying  the  energy  (all  matrix  elements 
are  energy  dependent)  until  a  zero  in  the  determinant  of  the  full  secular 
matrix  Is  found. 

The  most  severe  approximation  to  the  electronic  coupling  Interaction 
between  the  principal  and  environmental  subclusters  of  (1)  would  be  to 
assume  that  there  is  no  electronic  coupling  at  all,  i.e.,  that  Gab  and 
Gba  vanish.  This  would  put  (1)  in  block  diagonal  form  so  that  molecular 
orbitals  could  be  found  for  Ta_1  and  Tb"^  separately.  However,  these 


5 


orbitals  differ  from  those  found  by  calculations  of  the  isolated  sub¬ 
clusters  because  the  atomic  potentials  used  to  create  the  matrix  elements 
of  T  in  (1)  are  assumed  to  include  the  electrostatic  effects  of  the 

a 

atoms  in  the  environmental  subcluster.  We  call  this  level  of  approxima¬ 
tion  the  "molecular -field"  approximation. 

To  go  beyond  the  molecular -field  approximation,  the  effects  of 
electronic  coupling  between  the  subclusters  must  be  considered  as  well 
as  the  electrostatic  coupling.  The  method  presented  here  is  called 
"iterative"  coupling.  There  are  two  components  of  this  approximation: 
first,  contraction  of  the  secular  matrix;  second,  evaluation  of  the 
matrix  elements  of  the  environment  only  at  certain  trial  energies 
(to  be  discussed  further  below). 


Contraction  of  the  full  secular  matrix  Is  accomplished  using  a  simple 
property  of  linear  algebra  [5].  For  example,  the  linear  system  of  (1)  can 
be  re-written  In  a  way  that  eliminates  the  Afa  coefficients.  To  see  this, 
the  matrix  multiplication  of  (1)  Is  explicitly  written  out  to  obtain  the 
following  two  matrix  equations: 


A.  +  G,b 


V 


(2) 


Gba  A.  +  V1  Ab  *  0 


Solving  (3)  for  Ab  and  substituting  into  (2), 


ITa’  -  Gab  Tb  Gba)  Aa  3  0 


(3) 


(4) 
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where  Tb  is  the  matrix  inverse  of  T^.  Here  the  order  of  the  secular 
matrix  has  been  reduced,  i.e.,  the  matrix  has  been  "contracted";  solutions 
to  the  linear  system  of  (4)  are  obtained  only  when  the  determinant  of  the 
matrix  in  brackets  in  (4)  vanishes.  Once  the  column  vector  Aa  has  been 
obtained,  Ab  can  be  found  using  (3).  Thus,  solution  of  the  contracted 
linear  system  of  (4)  is  equivalent  to  solution  of  (1),  with  the  exception 
of  cases  in  which  the  matrix  inverse  of  Tb_1  does  not  exist  (i.e.,  when 
det  (Tb-1)  vanishes). 

Taking  the  determinant  of  the  smaller  matrix  of  (4)  is  certainly  more 
computationally  efficient  than  calculating  the  determinant  of  the  larger 
matrix  of  (1);  however,  (4)  was  obtained  by  taking  the  matrix  inverse  of 
Tb-1,  so  there  is  probably  no  overall  computational  advantage  in  following 
this  procedure.  However,  consider  a  state  which  is  primarily  localized  in 
the  principal  subcluster  (T”^);  for  such  a  state,  the  molecular- field 

a 

energy  should  be  close  to  the  corresponding  full  secular  matrix  energy. 
This  suggests  the  following  iterative  prescription  for  an  approximate 
coupling  method: 

(1)  First,  an  energy  eigenvalue  (EQ)  is  found  from 
submatrix  T  using  the  molecular-field  approximation. 

(2)  The  matrix  elements  of  Gab,  Tb_1,  and  Gfaa  are  calculated 
using  EQ  as  the  trial  energy;  the  matrix  inverse  of  Tb_1 
is  found. 


(3)  Determinants  of  the  matrix 

Y1(E)  '  Sab<Eo>VEo>  Gba<Eo>  ' 

are  calculated  at  various  values  of  E  until  a  zero 
is  found.  This  value  of  E  (call  it  E-j)  is  assumed 
to  be  an  approximate  solution  of  (1)  (or  (4)).  If 
E-j  differs  considerably  from  EQ,  the  process  is 
repeated  from  step  2,  using 

E  -  6E0  +  (1  -  S)  E-j  ( 

as  the  next  trial  energy  (0.0  <  8  <  1.0). 

Note  that  in  the  third  step  of  this  procedure,  the  matrix  elements  of 
Sab.  Tb"^*  and  Gba  are  unchan9ed,  so  they  must  be  calculated  only  once 
per  Iteration.  Similarly,  the  matrix  inverse  need  be  taken  only  once 
per  iteration.  Thus,  this  procedure  can  lead  to  very  substantial 
Improvements  In  computational  efficiency,  particularly  for  systems  in 
which  the  matrix  T^  is  large.  The  crucial  questions,  of  course,  are 
the  accuracy  and  speed  of  convergence  of  this  method.  These  questions 
will  be  Investigated  using  methanol  as  an  example. 

III.  Resul ts 

The  electronic  structure  of  methanol  can  be  easily  calculated  using 
the  partitioned  scattered-wave  method  without  contraction;  thus,  it  can 


serve  as  an  example  with  which  to  test  the  validity  of  the  iterative 
coupling  procedure.  Furthermore,  as  an  organic  molecule,  methanol  is  a 
prototype  for  molecules  which  might  need  the  iterative  coupling  procedure 
such  as  large  biological  molecules,  polymers,  or  embedded  clusters. 

The  partitioning  of  the  methanol  molecule  is  shown  schematically 
in  Fig.  1.  The  single  partitioning  sphere  containing  the  carbon  and 
oxygen  atoms  is  treated  as  the  principal  subcluster,  while  the  four 
hydrogen  atoms  collectively  form  the  environment.  The  starting 
potentials  for  the  iterative  calculation  were  obtained  from  a  fully 
converged,  partitioned  scattered-wave  calculation. 

The  second  column  of  energy  levels  in  Fig.  2  shows  the  results  of 
the  molecular-field  calculation  for  the  C-0  subcluster,  in  which  all 
scattering  from  the  four  hydrogen  atoms  is  neglected.  Note  that  two 
pairs  of  levels  are  degenerate:  the  la",  3a'  pair  and  the  2a",  5a' 
pair.  This  degeneracy  is  due  to  the  fact  that  the  true  symmetry  of  the 
CO  subcluster  is  C^,  so  that  the  degenerate  pairs  correspond  to  it  levels 
of  an  isolated  CO  molecule.  Even  though  the  potentials  in  this  molecular- 
field  calculation  include  the  electrostatic  effects  of  the  hydrogen  atoms, 
the  spherical  averaging  process  in  the  C  and  0  atomc  spheres  precludes 
any  symmetry  breaking. 

For  purposes  of  comparison,  a  standard  scattered-wave  calculation  of 
an  isolated  CO  molecule  was  performed  using  the  C-0  bond  distance  of 
methanol.  The  energy  levels  are  very  similar  to  those  of  the  molecular- 
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field  calculation,  except  that  they  are  shifted  lower  by  about  0.05-0.09 
rydbergs.  This  shift  is  explained  by  the  donation  of  charge  from  the 
hydrogen  to  the  carbon  and  oxygen  atoms  in  the  methanol  molecule. 

Because  the  molecular-field  calculation  uses  potentials  from  a  converged 
partitioned  scattered-wave  calculation,  the  energy  levels  reflect  the 
charge  transfer  but  not  the  stabilization  energy  associated  with  the 
interaction  with  the  hydrogen  atoms. 

The  third  column  of  Fig.  2  shows  the  energy  levels  of  the  f  ;t 
iteration  of  the  iterative  coupling  method.  The  initial  enviror  it 
energies  (EQ  in  (5))  were  taken  from  the  molecular- field  results 
Comparing  these  approximately  coupled  results  to  the  "exact"  or  fully 
coupled  results  in  the  fifth  column,  it  is  seen  that  although  the 
energy  levels  have  the  correct  ordering,  they  have  a  tendency  to 
"overshoot"  the  fully  coupled  values.  In  the  fourth  column  of  Fig.  2 
are  the  energy  levels  of  the  second  iteration;  these  used  initial 
environment  energies  obtained  from  the  average  of  the  molecular- field 
and  first  Iteration  energies,  corresponding  to  0  *  0.5  in  (6).  The 
energy  levels  resulting  from  the  second  iteration  are  in  remarkably 
good  agreement  with  those  of  the  fully  coupled  case. 

Plots  of  the  wavefunctions  of  the  3a'  level  from  molecular- field 
and  first  iteration  calculations  are  shown  in  Fig.  3  and  Fig.  4.  These 
can  be  compared  with  the  corresponding  wavefunction  of  the  fully  coupled 
calculation,  which  is  shown  in  Fig.  5.  The  wavefunction  from  the  second 
iteration  calculation  is  not  shown  because  it  is  virtually  identical  with 
that  from  the  fully  coupled  calculation. 
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IV.  Summary 

For  the  example  of  methanol,  the  iterative  coupling  method  successfully 
reproduces  the  energies  and  wavefunctions  found  in  the  fully  coupled  ("exact") 
calculation.  Yet,  the  iterative  coupling  method  involves  substantially  less 
computational  effort.  For  example,  the  fully  coupled  calcuation  of  the 
methanol  a'  state  uses  a  secular  matrix  of  order  93,  while  that  of  the 
approximate  calculation  has  order  68.  Because  the  numerical  calculation 
of  the  determinant  dominates  the  use  of  computer  time,  the  reduction  of 
the  order  of  the  secular  matrix  takes  on  increasing  importance  as  the 
cluster  size  increases.  (Calculations  involving  larger  clusters  have  been 
performed  and  will  be  reported  in  future  papers.) 

Although  fully  self-consistent  starting  potentials  were  used,  the 
iterative  coupling  calculations  reported  here  are  not  in  themselves 
self-consistent;  however,  there  is,  in  principle,  no  obstacle  to  ex¬ 
tending  the  method  in  order  to  perform  self-consistent  calculations. 

In  summary,  the  iterative  coupling  method  has  been  shown  to  be  an 
effective  way  to  drastically  reduce  the  computational  effort  involved 
in  scattered-wave  molecular-orbital  calculations.  Thus,  this  method 
has  the  potential  of  opening  up  new  areas  of  application  of  the 
scattered-wave  method,  applications  which  were  previously  impractical 
or  impossible  due  to  cluster  size.  These  include  biological  macro¬ 
molecules,  polymers,  and  microstructured  electronic  materials. 
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Fig.  1. 


Fig.  2. 


Fig.  3. 


Fig.  4. 


Figure  Captions 

Diagram  of  the  partitioning  of  the  methanol  molecule. 

The  single  partitioning  sphere  containing  the  carbon 
and  oxygen  atoms  is  treated  as  the  principal  subcluster, 
while  the  four  hydrogen  atoms  collectively  constitute 
the  environment. 

Molecular-orbital  energy  levels  from  iteratively  coupled 
scattered-wave  calculation  for  methanol  using  the 
subcluster  partitioning  of  Fig.  1,  compared  with  a 
standard  SCF-Xa-SW  calculation  for  an  isolated  CO 
molecule  on  the  left  and  the  "exact"  partitioning 
results  on  the  right.  The  second  column  shows  the 
molecular-field  energies  for  CO.  The  third  and 
fourth  columns  show  the  first  and  second  iteration 
approximate  coupling  results  for  methanol. 

Contour  map  of  the  3a'  (tt)  molecular-orbital  wavefunction 
for  CO  In  the  molecular- field  limit.  The  solid  and 
dashed  contours  represent  positive  and  negative  values, 
respectively,  of  the  wavefunction. 

Contour  map  of  the  3a'  molecular-orbital  wavefunction  for 
methanol  after  the  first  iteration  in  the  coupling  of  the 
hydrogen  environment  to  the  CO  subcluster.  Note  the 
formation  of  incipient  C-H  "a  bonds". 
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